Strongly correlated 2D quantum phases with cold polar molecules: 
controlling the shape of the interaction potential 
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We discuss techniques to tune and shape the long-range part of the interaction potentials in 
quantum gases of polar molecules by dressing rotational excitations with static and microwave fields. 
This provides a novel tool towards engineering strongly correlated quantum phases in combination 
with low dimensional trapping geometries. As an illustration, we discuss a 2D crystalline phase, 
and a superfluid-crystal quantum phase transition. 
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The outstanding feature of cold atomic and molecular 
quantum gases is the high control and tunability of micro- 
scopic system parameters via external fields. Prominent 
examples are the realization of low-dimensional trapping 
geometries, and the tuning of the contact interaction via 
Feshbach resonances Q . In this letter we extend this con- 
trol to the shape and the strength of interactions with 
the goal to generate new classes of potentials. In com- 
bination with low-dimensional trapping, these provide a 
framework for realizing new many body quantum phases 
and phase transitions. We elaborate these ideas in the 
context of polar molecules prepared in their electronic 
and vibrational ground states Q{. 

Cold polar molecules have attracted significant the- 
oretical interest, in particular in the context of dilute 
bosonic dipolar quantum gases 0, Q . Special focus was 
on the appearence of thermodynamic instabilities and ro- 
ton softening for weakly interacting gases [{|. Below we 
are specifically interested in polar molecules in the strong 
interaction limit, where the stability of the dipolar gas is 
guaranteed by a confinement of the particles into a two- 
dimensional (2D) setup. We show, that applying appro- 
priately chosen static and/or microwave fields allow us 
to design effective potentials V^P(R) between pairs of 
molecules (see Fig. da)). In turn, these potentials give 
rise to interesting new many-body phenomena. As an 
illustration, we consider the appearance of a crystalline 
phase, and an associated quantum melting to a superfiuid 
phase, for a 2D dipolar interaction VJP(R) = D/R 3 . 
The interaction strength is characterized by the ratio 
between the interaction energy and the kinetic energy 
rd = Dm/h 2 a with a the average interparticle distance. 
We determine the transition point tqm = 18 ± 4 of this 
superfiuid to solid quantum phase transition via Path 
Integral Monte Carlo simulations. We find that realistic 
experimental parameters of polar molecules allow for the 
realization of this quantum phase, which has never been 
observed so far with cold atomic or molecular gases. 

We consider heteronuclear molecules prepared in their 
electronic and vibrational ground states. We focus on 
bosonic molecules with a closed electronic shell = 




FIG. 1: (a) Effective potentials: dipole potential induced 
by a static electric field V^P = D/R 3 (solid line); poten- 
tial with a single microwave transition with Q/A = 0.1 
(dashed line); attractive potential induced by an additional 
microwave coupling to the states |1,±1) with Ai = 100A 
and fix /Ax = 0.3 (dash-dotted line), (b) Contour plot of 
the potential V(r) l\/D (Eq.Q) in the (R, 2)-plane. The in- 
stanton solution is shown as dashed line. 



0), e.g. of the type SrOjRbCs or LiCs, as most relevant 
in present experiments |2j . The Hamiltonian for a single 
molecule i is £T« = p 2 /2?« + E trap (r l ) + - d<E(t). 
The first two terms are the kinetic energy and the 
trapping potential for the center-of-mass motion of the 
molecule of mass m. The term H^ t describes the internal 
low energy excitations of the molecule, i.e., the rotational 
degree of freedom of the molecular axis. This term is well 



described by a rigid rotor H^ t = BJf with B the rota- 
tional constant (in the few to tens of GHz regime) and 
Ji the dimcnsionless angular momentum. The rotational 
states | J, M) with quantization axis z, and with eigenen- 
ergies BJ(J + 1) are coupled by a static or microwave 
field E via the electric dipole moment (typically of 
the order of a few Debye) . We assume a setup where the 
molecules are confined to a 2D configuration in the x, y- 
plane by a tight harmonic trapping potential midj_zf/2 
with u>± = fijma\, as provided by an optical poten- 
tial. We remark that the different dynamic polarizabili- 
tics parallel and perpendicular to the molecular axis give 
rise to tensorshifts, which provide an additional (small) 
state dependent potential for excited rotational states 0] . 
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The interaction of two polar molecules at distance 
r = |r| = |ri — T2| is described by the Hamiltonian 
H = J2i=i H + ^dd with the dipole-dipole interaction 
Vdd = [did 2 — 3(din)(d 2 n)]/r 3 and n = r/r as unit vec- 
tor . In the absence of external fields E, the interaction 
of two molecules in their rotational ground state is deter- 
mined by the van der Waals attraction V v dw ~ —Cg/r 6 
with Cg as d 4 /6B, valid outside of the molecular core 
region r > r rot = (d 2 /B) 1 ^. By dressing with an ex- 
tcrnal field E. we induce and design long-range interac- 
tions. On a formal level the derivation of the effective 
interaction proceeds in two steps: (i) We derive a set of 
Born-Oppcnhcimer (BO) potentials by first separating 
into center-of-mass and relative coordinates, and diago- 
nalizing the Hamiltonian for the relative motion for fixed 

molecular positions r. H re \ = J2i=i 2(^0! — d,E) + Vdd- 
Within an adiabatic approximation, the corresponding 
eigenvalues play the role of an effective V^P(r) interac- 
tion potential in a given state manifold dressed by the 
external field • (u) We eliminate the motional degrees 
of freedom in the tightly confined z direction to obtain 
an effective 2D dynamics with interaction VJP(R). In 
the following we consider the cases of a static field and a 
microwave field coupling the lowest rotor states. 

A static electric field applied perpendicular to the 
trap plane, E = E z e z , polarizes the rotational ground 
state of each molecule, and induces finite dipole moments 
(dj) = \/De z . These give rise to a long-range dipole- 
dipole interaction, V^P(r) = D(r 2 — 3z 2 )/r 5 , where D 
is tunable by the field strength E z . This expression of 
interacting dipoles aligned by the external field is valid 
for distances larger than r > r+ with r* defined by 
Ce/ r * ~ D/r 3 . Furthermore, for r > r rot the BO ap- 
proximation is easily fulfilled. The combination of the 
dipole-dipole interaction and the transverse trapping po- 
tential implies an interparticlc potential 
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with l± = (Dm/2H 2 a±) 1 / 5 a± (see Fig.^b)), and where 
we require lj_ max(r^, r rot ). The above potential 
has a saddle point at R s as 1.28Zj_ and z s as ±0.64Zj_ 
with the potential height V s as Q.34D/lj_. As a con- 
sequence, the short distance regime, r < l±, is sep- 
arated by a potential barrier from the large distance 
regime, r > Zj_. Within a semi-classical approxima- 
tion, the tunneling rate through this barrier takes the 
form T = ui p exp[— c(Dm/2h 2 a±) 2 ^ 5 ] with c as 5.86 and 
lo p ~ y^D/ma 5 the "attempt frequency" . Therefore, in 
the strongly interacting limit with a tight confinement 
the tunneling rate is exponentially suppressed, and the 
probability to find two particles at a distance r < l± is 
negligible, and guarantees the stability of an ensemble 
of polar molecules. Then, the system is completely de- 
termined by the interaction potential at large distances 



r ^> l±, which is in particular also much larger than the 
short distances scales r*, r rot . Thus we can reduce the in- 
teraction to an effective 2D potential by integrating out 
the fast transverse motion as 



V e 2 ff U (R) = J d Zl J dZ2W±(ziW\M**)rVg(T), (2) 

where ipJ-( z i) — cxp(—z 2 /2a 2 ± )/ (Tra 2 L ) 1,>4 is the ground 
state harmonic oscillator wave-function in the tightly 
confined z-direction, and is valid for distance R > 
(D/hw^) 1 / 3 (note, that l± ~ (D/hw±_) 1/3 for realistic 
parameters). Therefore, the effective 2D interaction re- 
duces to VJP ~ D/R 3 for large separations R^> l±. 

Microwave fields can drive the transition of a molecule 
from the ground to the first excited rotational state. We 
denote the corresponding Rabi frequency by Q and the 
detuning by A. In the weak driving limit O < A the 
molecules are essentially in their rotational ground state 
with AC Stark shifted energy and a small admixture 
from the excited state. For two molecules approach- 
ing each other and blue detuning A > (cf. Fig. 2), 
the microwave field will be resonant at distance Rq with 
d 2 /R,Q ~ HA, i.e., the ground and first excited rotational 
states will be strongly mixed. This gives rise to an ef- 
fective BO potential, which for r < Rq exhibits a strong 
repulsion inherited from the ~ 1 /r 3 excited state dipole- 
dipole interactions; Similar ideas have been discussed as 
"blue shield" in the context of alkali atoms Q- 
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FIG. 2: (a) Energy levels of the rotor in a weak electric 
field. The microwave transition with detuning A and Rabi 
frequency fl is shown as the arrow, (b) BO potentials for the 
internal states for Q, = 0, where \g\e)± = (|<?;e) ± \e;g))/y/2 
denote the gerade/ungerade first excited states, (c) Avoided 
level crossing due to the microwave coupling with the effective 
BO potential V e 3 ff D given by the highest eigenenergy state. 

As an example, we consider a microwave field with 
polarization along the z-axis, driving the transition be- 
tween the ground state and the first excited state with 
M = 0. Excitations into higher states are suppressed due 
to the anharmonicity of the spectrum, and spontaneous 
emission is irrelevant in the microwave regime. The in- 
teraction Hamiltonian also couples the excited state |1, 0) 
to the states |1,±1) for non-vanishing z (note that such 
a coupling can also be induced by the tensor shift). In 
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the following, we suppress this coupling by introducing a 
(weak) static electric field E z <C B / d lifting the degen- 
eracy of the J = 1 manifold by S ^ A, see Fig. |2Ia-b), 
and we focus on distances r > r$ = (d 2 /hS) 1 / 3 . For each 
polar molecule, the system reduces to a two level system 
with the ground state \g) and the first excited state |e) 
(M = 0) dressed by the electric field, and energy dif- 
ference hiu e . In a rotating frame the Hamiltonian H rc \ 
expressed in the basis {\g;g), \g;e), \e;g), |e;e)} becomes 



T < huj±, and obtain the many body Hamiltonian [Tj 

,2 
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where D g = \(g\di\g)\ 2 and D e = |(e|di|e)| 2 denote 
the weak dipole interaction induced by the static field 
Eg, while D m d 2 /3 accounts for the dipole transi- 
ton element between the ground and excited state, and 
v(r) = (r 2 — 3z 2 )/Hr 5 is the shape of the dipole inter- 
action. The effective BO potential V^P(r) derives now 
from the eigenenergy of the state adiabatically connect- 
ing to \g,g) for r — > oo, and is given by the highest en- 
ergy level, see Fig. |2Jc). We find that for a detuning 
A with R = (D/hA) 1/3 > l± > r s the combination 
of the transverse trapping potential with V^g provides 
again a large tunneling barrier separating the short dis- 
tance regime from the long distance regime. We em- 
phasize that the existence of this barrier makes the sys- 
tem stable - in contrast to the incomplete shielding dis- 
cussedpreviously in the context of "blue shielding" with 
lasers |jj . Then, in analogy to the discussion of the static 
electric field, we obtain the effective 2D potential via inte- 
grating out the transverse motion; the effective potential 
V 2 ^ is shown in Fig.n^a). The distance Ro separates a 
weakly interacting regime at large distances R > Rq with 
yjp ~ [(2fl 2 /A 2 )D + D g ] /R 3 from a strongly repulsive 
regime with V 2 ^ ~ D/R 3 on distances l± < R < Ro- 

The BO approximation above is valid as long as 
the passage through the level crossing at Ro is adi- 
abatic. For a realistic setup with average interpar- 
ticle distance a ~ Ro, the average particle velocity 
can be estimated as h/ Rom. Then, the Landau-Zcncr 
probability for a diabatic crossing reduces to P LZ = 
cxp [— 2ir(DmQ 2 / h 2 RoA 2 )] , and the BO approximation 
is valid for A 2 /f2 2 < Dm/h 2 Ro. This condition competes 
with the weak coupling constraint fl 2 /A 2 < 1. However, 
for f2/A ~ 0.1, both conditions can be satisfied for strong 
dipole interactions with > 100. 

The use of additional microwave fields coupling to the 
|1, ±l)-states with detuning A^ ^> A and Rabi fre- 
quency Qj_ allows for further shaping the effective po- 
tentials. One such example is shown in Fig. 1(a), where 
VJP in the long range part R > Ro becomes attractive 
and allows for the existence of bi-molecular bound states. 

Finally, we can extend the above discussion to a gas of 
cold polar molecules confined into 2D with temperature 
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The first term accounts for the kinetic energy within the 
x,7/-plane, while the second term denotes the effective 
interaction potential with the strong repulsion ~ D/R 3 
on short distances r < Ro, see above. The validity of 
this effective Hamiltonian requires that tunneling events 
through the barrier of the interparticle potential are sup- 
pressed, i.e., strong interactions and tight confinement 
with Dm/h 2 a± 3> 1. In turn, for decreasing interac- 
tions, tunneling through the barrier takes place and three 
body recombinations become relevant driving a crossover 
into a potentially unstable regime (note, that the ap- 
proach towards this intermediate region from the weakly 
interacting side has been previously discussed @). The 
Hamiltonian Eq. ||2J gives rise to novel quantum phenom- 
ena, which have not been accessed so far in the context of 
cold atoms/molecules. As an illustration, we focus on the 
interaction VjJ^R) = D/R 3 with the dimensionless pa- 
rameter Td = Dm/h 2 a and the particle density n = 1/a 2 . 
In the following, we present a tentative phase diagram of 
the system, see Fig. and show that for realistic ex- 
perimental parameters the polar molecules can be driven 
from the superfluid into the crystalline phase. 

In the limit of strong interactions = Dm/h 2 a ^> 1 
the polar molecules are in a crystalline phase for tem- 
peratures T < T m with T m w 0.09£>/a 3 H- The config- 
uration with minimal energy is a triangular lattice with 
spacing a L = (4/3) 1 / 4 a. The excitations are two lin- 
ear sound modes with characteristic phonon frequency 
ujp = y/ D I 'a b m. The static structure factor S diverges 
at a reciprocal lattice vector K, and S(K)/N acts as an 
order parameter for the crystalline phase. In the opposite 
limit of weak interactions ra < 1, the ground state is a 
superfluid (SF) with a finite (quasi) condensate. The SF 
is characterized by a superfluid fraction p s (T), which de- 
pends on temperature T, with p s (T = 0) = 1. For finite 
temperature, a Berezinskii-Kosterlitz-Thouless transi- 
tion towards a normal fluid is expected to occur at 
T KT = irp s h 2 n/2m. In turn, very little is known on the 
intermediate strongly interacting regime with > 1; 
see for a discussion in ID. Here, we focus on this 
intermediate regime and determine the critical interac- 
tion strength r QM for the quantum melting transition. 
We use a recently developed PIMC-code based on the 
Worm algorithm [lfj . In Fig. CSfd-e), the order param- 
eters p s and S(K.)/N are shown at a small temperature 
T = 0.014D/a 3 for different interaction strengths ra and 
particle numbers ./V = 36, 90. We find that p s exhibits a 
sudden drop to zero for as 15, while at the same po- 
sition S(K) strongly increases. In addition, we observe 
that in a few occasions p s suddenly jumped from to 
1, and then returned to 0, in the interval ~ 15 — 20, 
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FIG. 3: (a) Tentative phase diagram in the T — r<i plane: 
crystalline phase for interactions > t-qm and tempera- 
tures below the classical melting temperature T m (dashed 
line). The superfluid phase appears below the upper bound 
T < irh 2 n/2m (dotted line). The quantum melting tran- 
sition is studied at fixed temperature T = 0.014D/a 3 with 
interactions ra = 5 — 30 (dash-dotted line). The crossover to 
the unstable regime for small replusion and finite confinement 
uj± is indicated (hatched region) . (b) PIMC-snapshot of the 
mean particle positions in the crystalline phase for iV = 36 
at r c i ~ 26.5. (c) Density-density (angle-averaged) correlation 
function g2(r), for N = 36 at r& ~ 11.8. (d) Superfluid den- 
sity p s and (e) static structure factor S(K)/N as a function 
of rd, for N — 36 (circles) and N = 90 (squares). 

which suggests a competition between the superfluid and 
crystalline phases. These results indicate a superfluid to 
crystal phase transition at r QM = 18 ± 4. The step-like 
behavior of p s and S(K) /N is consistent with a first order 
phase transition. For > r QM the crystalline structure 
is triangular, in agreement with the discussion above, 
see Fig. Etb). Note, that the superfluid with ~ 1 
is strongly interacting, in particular, the density-density 
correlation function is quenched on lengths R < a, see 
Fig. Etc). Moreover, the (quasi) condensate involves a 
fraction of the total density, and therefore small coher- 
ence peaks in a time of flight experiment are expected. 
However, the detection of vortices can be used as a defini- 
tive signature of superfluidity, while Bragg scattering 
with optical light allows for probing the crystalline phase. 

While many of the molecular species studied at present 
in the lab are candidates for observing the above phenom- 
ena, we note that e.g. SrO has a large dipole moment 
d = 8.9D and allows optical trapping with a red de- 
tuned laser with wave length A ~ 1/rai. Then, for a trap- 
ping potential with a± ~ 50nm and interparticle distance 
a ~ 300nm, wc obtain ~ 410 3> r QM , and the tunnel- 
ing rate is suppressed with d 2 m/h 2 a± ~ 2500. Since the 



classical melting occurs at T m (rd = 400) ~ 2/iK, this 
molecule is a candidate to reach the quantum regime. 

In conclusion, we have shown that cold polar molecules 
allow to extend the control of interactions in cold gases 
to the shape of the inter-particle potential, which opens 
a route towards generation of (novel) strongly corre- 
lated quantum phases in low dimensional trapping ge- 
ometries. In particular, we have discussed the appear- 
ance and properties of a crystalline phase, and an asso- 
ciated quantum phase transition to a dipolar quantum 
gas. Besides the fundamental interest of observing crys- 
tal phases with neutral particles in atomic and molecu- 
lar physics, the properties of dipolar molecular crystals 
might lead to interesting applications: these crystals ap- 
pear in the high density limit, i.e. for small interparticle 
separation, which could potentially be smaller than op- 
tical lattice structures. At the same time, crystal phases 
are expected to be stable in the sense that (possibly dam- 
aging) short range collisions are suppressed. 
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